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ABSTRACT 

Context. The environment of high-mass X-ray binaries can be characterized either by the SNR that forms these systems, 
or by the wind from the companion massive star. These regions should be tenuous but very hot, and surrounded by a 
dense and cold shocked ISM shell. The interaction between the jet and such a complex medium, also affected by the 
system proper motion, can lead to very different jet termination structures. 

Aims. The evolution of the jet termination regions during the life of a high-mass microquasar is simulated to improve the 
present understanding of these structures. Also, the evolving emission characteristics are modeled to inform potential 
observational campaigns for this class of object. 

Methods. We have performed 2D numerical simulations of jets propagating in different scenarios, corresponding to 
different epochs after the formation of the high-mass X-ray binary, using the code Ratpenat. We have also made simple 
estimates of the non-thermal emission that could be produced in the jet termination regions. 

Results. We find that, in the way through the hot and tenuous medium of the shocked wind/SNR ejecta, the jet suffers 
recollimation shocks in which it loses part of its thrust and ends in a strong shock inflating a hot cocoon. The jet head 
propagates with a speed similar to the medium sound speed, until it eventually reaches the denser and colder shocked 
ISM and the unperturbed ISM later on. In these last stages of evolution, the jet is significantly slowed down and can 
be disrupted. For relatively old sources, the microquasar peculiar velocity becomes important, leading to complete jet 
destruction. Extended non-thermal radiation can be generated in the jet termination regions, being hard X-rays and 
TeV photons the best suited wavelengths to observe these structures. 

Key words. X-rays: binaries-ISM: jets and outflows-Radiation mechanisms: non-thermal 
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' 1. Introduction a strong shock will form at the tip of the jet, and particle 

, ... /TT^ /r^r-r^x 1 1 acceleration may tap a significant fraction of the jet lumi- 

1^ High-mass X-ray bmaries (HMXB) are the parent popula- ^osity there 

. tion of high-mass microquasars (HMMQ) , the latter having ^^^^ ^^^U^^ ^.^^ dynamically negligible, the jet 

^ the pecuharityofpresenting radio jet actmty (e.g. Mirabel propagate as far and fast as the inertia of the swept ISM 

O & RodriguesHMl Ribo 2005). Jets o HMMQs propagate p^.^i^s and, eventually, the jet head is shocked inflating a 

^ m a medium characterized by the stellar wind of the mas- ^.^^ ^^^^^^ ^^^^ j^^^ ^ ^j^^U ^^^^^^^ jg^ ^j^^^^^ 

sive primary star of the system At the scales of the binary, 2009), as in FRII galaxies (e.g. Fanaroff & Riley 1974; 

^> the stellar wind can significantly affect the dynamics of the j^^.^^^ ^ Alexander 1997). However, for systems hosting 

^ jet, trigger shocks suitable for particle acceleration, and un- Qg ^^^^^^ dynamics of the jet may be severely affected 

. der some conditions, disrupt the jet by a combination of ^^^^^^ ^^^^j^.^ unperturbed ISM. 

^ • asymmetric recollimation shocks and hydrodynamical m- ^^^^^ ^^^^^^ understand the hydrodynamical jet 

• - - stabilities see Perucho & Bosch-Ramon 2008 and Perucho ^^^^^^^^^ different stages of a HMXB, we have performed 

^* , /* , , „ , , . , 2D simulations of HMMQ jets for three different scenar- 
When launched, powerful enough jets can traverse the ^^^^ ^^^^^ ^j^^ g^^j^ explosion (case 1), the phase of 
companion wind, and once at distances from the injection ^j^^d/ISM interaction (case 2), and the phase in which the 
point much arger than the separation distance (z » dorb, ^ind/ISM interaction is affected by the system proper mo- 
where ^ is the axial coordinate and distance to the com- ^^^^^ j^^^^ ^1^^ estimated the broadband non- 
pact object), they do not suffer significant interference ^^^^^^^ emission that may be produced in cases 1 and 2, 
by the supersonic wind, which comes from behind the and made predictions of their main observational features, 
jet. At this stage, the jet head propagates with a speed 
^ nj ((iorb ) / (<^orb ) '^j ~ 0-1 '^j assumiug free expansion, 

where f j is the initial jet velocity. In this propagation phase, 2. The physical picture 

Send offprint requests to: V. Bosch-Ramon, e-mail: For a young HMXB, say of age t^rc ~ 3 X 10"^ yr, with typ- 

vbosch@mpi-hd.mpg.de ical ISM densities njsM ~ 1 cm~^ and moderate peculiar 
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velocities < 10^ cm s~^, the system may still be em- 
bedded inside an adiabatically expanding SNR (case 1). In 
such a situation, if a jet is launched, after ~ 10 yr it will 
encounter a region of shocked stellar wind confined by the 
shocked SNR ejecta, the latter to its turn driving a strong 
shock in the ISM. Once in the high-pressure region of the 
shocked stellar wind, the jet head will be slowed down to 
medium subsonic velocities (~ 1000 km s~^ for typical jet 
and wind properties, see below). At farther distances, the 
jet head will encounter the shocked SNR ejecta eventually 
reaching the shocked ISM. Finally, the jet will enter the 
normal ISM having lost significant amounts of thrust, pro- 
ducing a shock in the ISM and possibly suffering disrup- 
tion. Intermittent jet activity with low duty cycle will not 
clean the way for later ejections, since the sound speed in 
the shocked wind and the SNR ejecta is high. High enough 
duty cycles 10%) for the jet to keep clean its way could 
be understood in terms of an averaged persistent activity. 
The case of a HMMQ persistent jet interacting with a SNR 
has a clear example in SS 433 (see Velazquez & Raga 2000, 
Zavala et al. 2008). 

Older HMXBs have already left the SNR, and the stellar 
winds can directly interact with the ISM (case 2). In this 
case, the stellar wind inflates a bubble filled with shocked 
wind material, and this hot bubble drives a shock in the 
ISM that becomes radiative even for moderate ages of the 
source {tsrc ^ 10^ yr; see Castor et al. 1975). Therefore, a 
jet launched in this HMXB phase will interact first with the 
shocked wind bubble, reaching later on a thin, dense and 
cold shell of shocked ISM. Qualitatively, the situation is 
similar to the previous case, but now the ISM shell is much 
denser because of radiative cooling. Eventually, the jet can 
reach the unperturbed ISM, shocking it. As in case 1, the jet 
may be eventually disrupted. Cygnus X-1 could correspond 
to the situation represented by case 2 (Marti et al. 1996; 
Gallo et al. 2005, Russell et al. 2007). Since there are not 
clear signatures of a reverse shock in this source (Bordas 
et al. 2009), and only the shocked ISM thermal radiation is 
detected, it may be that the Cygnus X-1 jet head is fully 
disrupted. A shock driven in the ISM can still develop in 
such a conditions since the disrupted jet region has much 
higher pressure than the environment as simulations of FRI 
radio galaxies show (Perucho & Marti 2007; see also Bordas 
et al. 2010a). 

If the HMXB age is old enough, say t^rc ^ 10^ yr, 
the wind-driven ISM shock can be slower than v^. Then, 
a bow-shaped shock forms, driven by the system proper 
motion and pointing in its direction (case 3). In this sce- 
nario, the shocked ISM material is slowly evacuated in the 
opposite direction to the system motion, and the shocked 
wind material is also advected in the same direction but 
at much faster velocities, similar to shocked wind sound 
speed (~ 10^ km s~^). If the jet is launched in this stage, 
for jet directions misaligned enough with the binary mo- 
tion the shocked wind material will impact the jet roughly 
from the side. At farther distances, the jet will encounter 
the shocked ISM, roughly as dense and cold as in case 2. 
The dense shocked ISM layer plus the shocked wind lat- 
eral impact will stop the jet advance and lead to total jet 
destruction before reaching the unshocked ISM. 

Sketches of the scenarios corresponding to cases 1, 2 
and 3 are presented in the top, middle and bottom panels 

of Fig. m 
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Fig. 1. Sketches of a jet propagating through: the shocked 
stellar wind, the shocked SNR ejecta, and the shocked 
and unshocked ISM (top); the shocked stellar wind, and 
the shocked and unshocked ISM (middle); through the un- 
shocked and shocked stellar wind, and interacting with the 
shocked ISM under a significant peculiar velocity (bottom). 



3. Studying the jet termination 

3.1. Medium characteristics 

We have adopted an analytical approach to characterize 
the medium in which the jet propagates. In the three cases 
considered here, the ISM density has been fixed to njsM = 
1 cm~^, with a temperature Tjsm = 10^ K, as expected in 
regions close to massive stars. The stellar wind velocity and 
the stellar mass-loss rate have been taken to be 'L'w = 2 x 
10^ cm s~^ and = 10~^ Mq/jt, respectively (typical for 
an O star; see, e.g.. Puis et al. 2009), which corresponds to a 
wind kinetic luminosity = M^v1^/2 ^ 1.3x10^^ ergs~^. 

In case 1, the system is assumed to be at the cen- 
ter of a SNR. The ejected mass of the SNR has been 
taken to be 2 M©, and the explosion energy 5 x 10^^ erg. 
These moderate values have been adopted accounting for 
the fact that the energetics of supernovae leading to the 
formation of some HMXBs could be quite low (a rela- 
tively extreme case may be Cygnus X-1; see Mirabel & 
Rodrigues 2003). We have calculated the properties of the 
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shocked ejecta of the SNR and the shocked ISM adopt- 
ing the solution of a spherical explosion in its adiabatic 
phase (e.g. Blandford & McKee 1976) for tsrc 3 x 10"^ yr. 
The radius of the whole region affected by the SNR is 

^bs ~ {EsNR/munisuy^^ts/c ~ 3 x 10^^ cm for the 
adopted parameter values, where mn is the hydrogen mass. 
The density of the shocked ISM is the one for a strong 
adiabatic shock, njsMs = 4nisM5 and the temperature, 
^iSMs ^ 4 X 10^ K. From Zbs and ^snr, we derive a charac- 
teristic pressure within that radius, Psnr ~ 2 £^snr/3 Vbs ~ 
3 X 10-^0 erg cm-^ (Tsnr ^ 10^ K; ^bs = 4 tt z^JS), as well 
as the locations of the contact discontinuity between the 
shocked ISM and the shocked ejecta, ~ 0.9 Zbs- The loca- 
tion of the contact discontinuity between the shocked stellar 
wind and the SNR ejecta is computed from taking the same 
pressure in both sides, obtaining Zws ^ 8 x 10^^ cm. The 
shocked wind density is n^s ^0.02 m~^, and the temper- 
ature Tws ^ 10^ K. We simulate the jet once it is already 
inside the shocked wind region, since at closer distances the 
dynamical impact of the medium on the jet is negligible. 

In case 2, we have adopted the solution given by Castor 
et al. (1975) that describes the interaction of a stellar wind 
with the ISM at the stage when the shock in the latter 
is already radiative. We have adopted in this case an age 
of 10^ yr, and the radius of the whole interaction region 
is Zbs = 0.76(Lw/mHnisM)^/^t3/^ ^ 2 x 10^^ cm. The 
pressure in the shocked wind region is Pws ^ ^ISM i^bs ^ 
3 X 10-^^ cm s-^ (Tws « 6 X 10^ K), and the density, 
nws = 4 X 10~^ cm~^. The shocked ISM region is a very 
thin, cool and dense layer. For simplicity, the width of this 
region has been taken equal to one simulation cell: 10^^ cm, 
with nisMs ~ 10^ cm~^ and Tjsms = 2 x 10^ K. As in case 
1, the jet injection point is already in the wind shocked 
region. 

In case 3, taking = 10^ cm s~^ and tgrc = 10^ yr, 
one gets ibs ^ lO^cms"^ <^ Vp and the wind-driven 
shock in the ISM becomes strongly asymmetric, with the 
shocked wind material advected backwards with respect to 
this shock. The wind driven shock in the ISM is now closer 
to the HMMQ than in case 2, at ^ 8 x 10^^ cm, at 
which the (pre-shock) stellar wind density and tempera- 
ture are = 2.4 x lO""^ cm"^ and = 10^ K. The 
density, pressure and temperature in the shocked wind 
are nws = 10"^ cm"^, Pws = 1-4 x 10"^^ erg cm"^ and 
Tws = 10^ K, respectively, and a lateral shocked wind ve- 
locity Vws± = 10^ cm s~^. The shocked wind region reaches 
Zws = 10^^ cm, in which a thin shocked ISM shell, of 
10^^ cm width, njsMs = 100 cm~^ and Tjsms = 10^ K, 
is located. In this simulation, the jet is launched in the un- 
shocked wind region, because it is more extended than in 
previous cases, and the unshocked wind is taken to come 
from behind the jet. 

The values of the main ambient medium parameters, for 
the three cases considered here, are summarized in Table 1. 

3.2. Simulations 

We have performed 2D simulations of jets evolving in the 
three different scenarios described above. In cases 1 and 2 
the jets are simulated using an axisymmetric grid in cylin- 
drical coordinates, whereas the transversal motion of the 
shocked wind in case 3 requires a slab jet in planar coordi- 
nates. We expect from our previous experience (see Perucho 



Case 1 Jet height medium dens. medium temp. 

z (cm) n {nip cm~^) T (K) 

Shocked wind 3 x 10'^ 002 10^ 

SNR 8 x 10^^ 0.02 10^ 

Shocked ISM 3 x 10^^ 4.0 4 x 10^ 

ISM 3.3 X 10^^ LO 10^ 

Case 2 

Shocked wind 3 x lO'"" 4 x 10"^ 6 x 10' 

Shocked ISM 2 x 10^^ 100 2 x 10^ 

ISM 2.01 X 10^^ LO 10"^ 

Case 3 

Stellar wind 1.25 x lO'"" 2.4 x 10"^ 10^ 

Shocked wind 8 x 10^^ 10"^ 10^ 

Shocked ISM 1 X 10^^ 100 10^ 

ISM 1.01 X 10^^ 1.0 10^ 



Table 1. The ambient medium in cases 1, 2 and 3. 



& Bosch-Ramon 2008, Perucho et al. 2010) that the results 
of the slab simulation will be qualitatively correct. The jets 
are injected with power Lj = 3 x 10^^ erg s~^ at zq = 
3 X 10^^ cm for cases 1 and 2, and zq = 1.25 x 10^^ cm for 
case 3. In case 3, the power corresponds to a jet square sec- 
tion (i.e. the power crossing a surface equal to the jet width 
square). In the cases 1 and 2, the jet is assumed to have a 
radius at the injection point R^o = 0.1 zq. In case 3, the jet 
half-side length has been taken as 1.5 x 10^^ cm. We have 
fixed the jet velocity at = 10^^ cm s~^, the Mach number, 
Mj ^ 18, and the jet temperature, Tj = 1.2 x 10^ K at in- 
jection in all three cases. The jets are composed by protons 
and electrons and have densities of pj = 1.8 x 10~^^ g cm~^ 
for cases 1 and 2, and pj = 8.3 x 10~^^ g cm~^ for case 3. 
The adiabatic exponent for the jets, with the given temper- 
ature and composition is F ^ 1.58. The ambient medium 
in the different simulations follows the description given in 
Sect. 13.11 We note that the shock driven by the jet head in 
the medium (wind, SNR and ISM) is called here bow shock 
(which is more a sound wave rather than a real shock in 
the hot shocked wind and SNR ejecta), and the jet shocked 
region is called cocoon, following Bordas et al. (2009). 

The simulations have been carried out with the finite- 
volume code Ratpenat^ which solves the equations of rel- 
ativistic hydrodynamics in conservation form using high- 
resolution-shock-capturing methods. Ratpenat was paral- 
lelized with a hybrid scheme with both parallel processes 
(MPI) and parallel threads (OpenMP) inside each process 
(see Perucho et al. 2010). The simulations were performed 
using the computational facility Tirant^ at the ''Servei 
dlnformatica^^ of the ^'Universitat de Valencia^\ with 12, 
8 and 4 processors for cases 1, 2 and 3, respectively. In 
the axisymmetric simulations, the numerical grid box ex- 
pands transver sally to 100 i?jo, with an extended grid of 
100 i?jo added to the uniform grid. In the axial direction, 
the grid extends to 150 i?jo in case 1, and to 100 i?jo in case 
2. The numerical resolution of these simulation is of eight 
cells per initial jet radius. The computational grid, includ- 
ing the extended region is 1200 x 1000 cells. In the case of 
the simulation of the slab jet, the grid expands 100 Rjo to 
each side of the axis, plus 100 i?jo extra at each side, and 
100 i^jo in the axial direction. The resolution used here is of 
four cells per jet initial half-side length. The dimensions of 
the grid are thus 1200 x 400 cells. The boundary conditions 
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adopted in the Cartesian and cyHndrical cases are reflection 
at the jet bottom side, injection right at the jet base, and 
outflow in the rest. Reflection in the jet axis has been also 
adopted in the cylindrical simulations. 

3.3. Results 
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Fig. 2. Axial cuts of density (top, left), pressure (top, 
right), temperature (down, left) and velocity (down, right) 
at the last snapshot of the evolution in case 1 (tgrc = 
8900 yr). Units are cgs. 



The simulation of case 1 has been run during 8900 yr 
of evolution. During the first 4500 yr, the jet propagates 
up to z ~ 3 X 10^^ cm, with an advance velocity ^ 2 x 
lO^cms"^. The jet is initially underpressured with respect 
to the hot ambient medium (shocked wind). This generates 
a strong conical recollimation shock in the jet soon after 
injection, followed by pinching. At z ~ 3 x 10^^ cm, the jet 
head reaches the shocked ISM and its advance is stopped. 
The speed of the head is reduced to 2 x lO'^cms"^ and 
the jet is only able to propagate up to z ^ 3.2 x 10^^ cm 
by the end of the simulation. The expansion of the bow 
shock during this phase is mainly radial (perpendicular to 
z)^ with an expansion to a distance of 4.3 x 10^^ cm from 
the jet axis (outside the frame of the image). The reverse 
shock becomes stronger after the collision with the shocked 
ISM, and pinching is enhanced along the jet. By the end 
of the simulation, the SNR shock in the ISM has expanded 
farther, relaxing the density and pressure conditions that 
the jet finds at arrival, which makes easier for the jet to 
carve its way across this region. At the given velocity, the 
jet would need about 3000 yr more to complete its path 
through the shocked ISM. 

In case 1, the pressure in the cocoon and bow-shock re- 
gions is basically constant from the early stages of the sim- 
ulation: ^ 5 X 10~^^ erg/cm^, whereas the density starts 
increasing when the jet head reaches the ISM shocked by 
the SNR, from 3 x 10"^^g/cm^ to 5 x 10"^^g/cm^ during 
the period 4500-8900 yr. The temperature in the shocked 
ambient (shocked wind and SNR ejecta) falls slowly during 
the whole run from 2 x 10^ to 1.2 x 10^ K. In the cocoon, 
the pressure is very similar to that in this shocked ambi- 
ent. The cocoon density increases from 3 x 10~^^g/cm^ 
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Fig. 3. From top to bottom, maps of density, pressure, tem- 
perature and axial velocity at the last snapshot of the evo- 
lution in case 1 (tgrc = 8900 yr). Units are cgs. 

to 7 X 10"^^g/cm^ during the first 2000 yr due to the 
jet low pressure with respect to the ambient, and remains 
fairly constant afterwards due to the recollimation and 
pinching shocks. The temperature in the cocoon decreases 
from 3 X 10^^ to 10^^ K during these first 2000 yr, and 
becomes constant from that moment on. Turbulence and 
mixing of jet and external material takes place in the co- 
coon. Turbulence follows, down to the scales allowed by the 
simulation resolution, a spectrum similar to Kolmogorov, 
and the entrained external material actually dominates the 
mass content of the cocoon. The axial cuts of density, pres- 
sure, temperature and velocity at the last frame, presented 
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Fig. 4. Evolution of pressure (top), density (middle), and 
temperature (bottom) in the bow shock (left) and cocoon 
(right) regions for case 1. Units are cgs. 
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Fig. 5. Position of the jet head (left) and mean radius of 
the bow-shock region (right) for case 1. Note that the bow 
shock has already left the shown grid by the end of the 
simultion. Units are cgs. 
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Fig. 7. Same as in Fig. [3] but for case 2 (tgrc = 9800 yr). 
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Fig. 6. Same as in Fig. [2] but for case 2 (tgrc = 9800 yr). 



in Fig. [21 show the rich structure in the jet after the jet is 
pinched by the ambient pressure, with plenty of shocks, up 
to the point where the jet material reaches the strong re- 



verse shock in the interaction with the ISM shocked by the 
SNR. The recollimation and pinching shocks appear early 
in the simulation given the low initial lateral pressure of 
the jet with respect to the medium, the pinching structures 
getting enhanced when the jet head reaches the shocked 
ISM. The maps for the density, pressure, temperature and 
axial velocity at the last snapshot (8900 yr) of the evolution 
are shown in Fig. [3l The temporal evolution of the average 
cocoon and bow-shock density, pressure, temperature, and 
velocity, are shown in Fig. SI The temporal evolution of the 
jet head position and the mean radius of the bow shock are 
shown in Fig. [5l 

The simulation of case 2 has been run during 9800 yr 
of evolution. The jet crosses the stellar wind and shocked 
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Fig. 8. Same as in Fig. S] but for case 2. 




Fig. 9. Same as in Fig. [5] but for case 2. 



wind region, from injection at z ~ 3 x 10^^ cm to the con- 
tact discontinuity with the shocked ISM at z ~ 2 x 10^^ cm, 
at 1500 yr after injection, at a mean velocity of ^'h ^ 
3.5 X lO^cms"^. As in case 1, the jet is initially under- 
pressured and collimated by the ambient, which produces 
pinching oscillations along the jet body. After reaching the 
shocked ISM, the jet head gets flattened due to the large 
density wall. At that point, the pressure in the bow-shock 
and cocoon regions of the jet starts to decrease slower with 
time. The same happens with the density in the cocoon, 
whereas the density of the shocked ambient starts to in- 
crease, since the head of the jet interacts with a denser 
medium. The backflow expands sideways, fllling the sim- 
ulated region, whereas the head of the jet acts as a drill 
on the shocked ISM region due to its much higher pres- 
sure and crosses it at t 4300 yr after injection. Thus, it 
needs 3000 yr to cross this small shell. As in case 1, the 
internal structure of the jet becomes very rich due to the 
interaction between the previous pinching patterns and the 
backward moving waves produced at the shock of the jet 
flow with the stagnation point, at z ~ 2.4 x 10^^ cm. The 
axial cuts of the jet, presented in Fig. [6l show this effect 
and the increase of the temperature and decrease of the 
velocity after each collimation shock, in which strong pres- 
sure jumps are produced. After the jet head crosses the 
shocked ISM layer, it propagates in the ISM with a super- 
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Fig. 10. Maps of density (top) and pressure (bottom) at the 
last snapshot of the evolution for case 3 {t^rc ~ 2000 yr). 
The shocked wind moves from top to bottom. 



for the density, pressure, temperature and axial velocity at 
the last snapshot (9800 yr) of the evolution are shown in 
Fig. [71 The temporal evolution of the average cocoon and 
bow-shock density, pressure, temperature, and velocity, are 
shown in Fig. [H The temporal evolution of the jet head 
position and the mean radius of the bow shock are shown 
in Fig. [9l Overall, the jet evolution of case 2 is similar to 
that of case 1, with the cocoon also affected by turbulence 
and strong mixing. 

Simulation 3 has been run during ^ 2000 yr. The jet is 

18 , 



sonic speed of ~ 2.4 x 10 cms , an order of magnitude injected in the grid in the wind region, at z = 1.25x10 cm. 
smaller than that in the shocked wind region. The maps When the jet reaches the shocked wind region, at z = 
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Fig. 11. Maps of temperature (top) and axial velocity 
(bottom) at the last snapshot of the evolution for case 3 
(tsrc ^ 2000 yr). The shocked wind moves from top to bot- 
tom. 



8 X 10^^ cm, it gets the impact from the shocked wind that 
has been deflected by the propagation of the binary system 
in the ISM (see Fig.[TJ bottom). This wind is enough, in the 
slab simulation, to distort the bow shock of the jet and dis- 
rupt the latter as follows. After crossing the wind/shocked 
wind discontinuity, the bow shock is deviated to the right 
from injection in the images due to the shocked wind com- 
ing from top to bottom. The backflow is also deviated in 
this direction due to the generated pressure gradient. This 
generates an excess of backflow material in the unshocked 
wind region at the right of the jet, resulting in a devia- 
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Fig. 12. Same as in Fig. [10] but for a source age of 300 yr. 



tion of the jet to the left. The combination of this process 
with the shocks generated at the different discontinuities 
destroy the jet even inside the wind region. When the head 
of the jet reaches after 300 yr the shocked ISM region, at 
z = 10^^ cm, the jet has not enough thrust to penetrate 
this layer and all the injected flow is thus converted into a 
backflow that fllls the inner regions and also the unshocked 
wind region. The combination of disruption and advection 
does not allow the jet to reach the normal ISM regardless 
the time elapsed, and thus jets with the given power will 
be trapped inside the wind/shocked wind region, and ad- 
vected backwards with respect to the system motion. The 
maps for the density, pressure, temperature, and axial ve- 
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Fig. 13. Same as in Fig. [TT]but for a source age of 300 yr. 

locity, at the last snapshot (2000 yr) and after 300 yr of 
evolution, are shown in Figs. [TOl [TTJ [12] and [131 

4. Non-thermal activity 

In order to estimate what may be the potential non-thermal 
luminosity of the radiation produced when the jet head 
propagates within the SNR or the shocked stellar wind, we 
have used a model for the cocoon emission similar to that 
presented in Bordas et al. (2009). In short, the model as- 
sumes the cocoon as an homogeneous emitter, which is true 
to first order provided that the flow is roughly sonic. Our 
simulations show turbulent motion and jet/medium mixing 
in the cocoon, and the flow becomes trans-sonic far from 



the jet tip, but at this stage we will neglect this effect. In 
the radiation calculations, following Bordas et al. (2009), 
the magnetic field {B) energy density iu^) has been as- 
sumed to be a 10% of the matter energy density, and the 
non-thermal injection luminosity, a 1% of the jet power. 
Postshock magnetic field values close to, or not far below, 
equipartition are expected if the jet keeps a non- negligible 
fraction of its total energy as magnetic energy, which is 
a reasonable assumption provided that jets are magneti- 
cally launched (e.g. Komissarov et al. 2007; Spruit 2010). 
Note that for an optically thin emitter, the radiation lumi- 
nosity will be just proportional to this fraction. Unlike in 
Bordas et al. (2009), here only the cocoon has been con- 
sidered as a powerful emitter, since the recollimation and 
reverse shocks suffered by the jet are strong, but the bow 
shock is too weak to accelerate particles. The results for the 
cocoon radiation are presented in Sect. 14. ll in which we fo- 
cus on electrons only given the low cocoon densities that 
prevent efficient gamma-ray production from proton-proton 
interactions. These low densities imply also that relativistic 
Bremsstrahlung will be inefficient there. 

When the jet head impacts the shocked ISM thin shell, 
the highest energy particles accelerated in the cocoon re- 
gion may reach this shell, which in case 2 reaches densities 
^iSMs ^ 100 cm~^. Given the relativistic Bremsstrahlung 
(electrons) and proton-proton interaction (if protons are 
also accelerated) timescales and the source age, about a 1- 
10% of the luminosity in relativistic particles reaching that 
shell could be radiated at very high energies. We note that 
denser media would increase the efficiency, which is simply 
oc nisMs- 

For the case when the jet has reached the undisturbed 
ISM, if it is not disrupted, the non-thermal behavior can 
be described as in Bordas et al. (2009). This would apply 
to the emission produced during the second half of case 2 
simulation. We have not taken into account in this work 
particle acceleration in the turbulent and sheared flow in 
the recollimated jet and cocoon regions, but stochastic and 
shear acceleration may be relevant there (see Rieger et al. 
2007 and references therein). Some thermal radiation may 
be produced in the shocked ISM. This thermal component 
tends to be quite extended and, given the shock velocities, 
would peak in the UV. 

We do not consider here case 3 since, due to jet full 
disruption, conditions seem less suitable for particle accel- 
eration. It might be that stochastic and shear acceleration 
may still take place in the jet disrupted medium, although 
estimating the non-thermal emission for this case requires 
an in-depth study that is beyond the scope of this work. 



4.1. Cocoon emission: cases 1 and 2 

The simulation results are used to characterize the phys- 
ical properties of the shocked regions. The position and 
velocity of the bow shock, the internal energy density (iith; 
and hence ub = O.li^th), the temperature, and the rest 
mass density are used at any source age to account for the 
evolution of the particle energy distribution, which is af- 
fected by both adiabatic and radiative losses: relativistic 
Bremsstrahlung (negligible), inverse Compton (minor) and 
synchrotron (dominant). For the IC emission, the radiation 
field provided by the companion star dominates over that of 
the CMB. The star is taken to be a bright O, with luminos- 
ity 10^^ erg s~^ and temperature 4 x 10^ K. At a distance 
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Fig. 14. Synchrotron and IC SEDs for both simulated cases 1 (black/thick hnes) and 2 (blue/thin hnes) computed for 
tsTc = 10^ (left), tsrc = 3 X 10^ (middle) and tsrc = 9 x 10^ yr (right). 



Table 2. Model parameters for the shell and cocoon for 
three different source ages. 



Parameter value 
Jet power [erg s"'] S^TlO^ 
Non-thermal power fraction 0.01 

stellar field energy density {case 1) [erg cm~^] 3.1 x 10""^^ 

stellar field energy density {case 2) [erg cm~^] 5.2 x 10""*^^ 
magnetic field {case 1) [G] 4.1 x 10~^ 

magnetic field {case 2) [G] 1.4 x 10~^ 

Acceleration efficiency 0.02 
Injection power- law index 2.1 
Maximum energy {case 1) 10"^ TeV 

Maximum energy {case 2) 1.8 x 10^ TeV 



of ~ 10 pc, the energy density ratio is u^/uqub ^ 7, in- 
creasing towards the central engine like ex 

We have assumed that the recollimation and reverse 
shocks accelerate particles, simplifying them as just one 
accelerator with average properties. The shock velocity has 
been taken '^g '^j, which yields an acceleration efficiency 
E - (l/27r) {vJcfqBc ^ 0.02 x qBc (e.g. Drury 1983). 
As mentioned before, the assumed non- thermal injection 
luminosity is a 1% that of the jet, i.e. 3 x 10^^ erg s~^, and 
the magnetic to internal energy density ratio is 0.1. This 
gives a 5-strength of about ~ 4 — 6 x 10~^ G depending 
on tsrc These 5- values are one order of magnitude higher 
than those inferred by Tudose et al. (2006) and Safi-Harb & 
Petre (1999) for the jet termination regions in Circinus X-1 
(a low-mass microquasar) and SS 433, respectively. On the 
other hand, the ages of these sources are also about one 
order of magnitude older than those considered here, and 
thus a lower magnetic field is naturally expected in the for- 
mer cases. Under the adopted magnetic field, synchrotron 
losses are dominant in the cases studied here. For the condi- 
tions assumed, maximum energies are typically - 10^ TeV. 
We remark that the adopted non-thermal fraction is quite 
conservative, and there is room to increase it. The values 
of the radiation model parameters are listed in Table [2] for 
cases 1 and 2 when the source age is ~ 10^ yr. 

The spectral energy distributions (SED) of the non ther- 
mal synchrotron and IC emission for cases 1 and 2, for evo- 
lution times tsrc = 10^,3 x 10^ and 9 x 10^ yr, are shown 



in Fig. [m Depending on the age of the source and the 
case, and given that ub ^ i^*, the synchrotron luminosity 
can be between ~ 10 — 1000 times higher than the IC one. 
This difference is more pronounced in case 1. This is a re- 
sult of the smaller amount of energy pumped by the stellar 
wind compared to that of a SNR, which leads to a lower 
pressure and therefore a lower B in case 2. We note that 
the IC dominance may be only possible if the 5-strength 
were well below equipartition. At the end of the simulations 
(^src = 9 X 10^ yr), the synchrotron fluxes are, for a 3 kpc 
distance source, of ~ 20 mJy (case 2) and 100 mJy (case 

1) , and ~ 4 (case 2) and 7 x 10""^^ erg cm~^ (case 1), 
in radio and X-rays, respectively. The IC emission yields 
fluxes of ~ 6 X 10~^^ (case 1) and 10~^^ erg cm~^ s~^ (case 

2) , and ~ 5 x 10"-^^ (case 1) and 3 x 10"-^^ erg cm~^ 
(case 2), in the ranges 0.1 - 100 GeV and > 100 GeV, re- 
spectively. We recall that these fluxes might be easily one 
order of magnitude higher, since the non- thermal fraction 
is not much constrained. 

It is worth noting that the cocoon is in pressure equi- 
librium with the medium. Since the dynamical timescales 
of the SNR/ISM and wind/ISM interaction structures are 
longer than those generated by the jet, the magnetic field is 
likely to be approximately constant with time. In this con- 
text, the radio emission grows with the source age due to 
particle accumulation but the X-ray fluxes do not change 
much, because X-ray synchrotron emitting electrons have 
time to cool down from quite early times. Concerning IC, 
the increase of the available cooling time due to longer tgrc 
compensates the dilution of the stellar photon field (and 
subsequent reduction of the emission efficiency) produced 
by the increase in cocoon size. This renders fairly constant 
total IC fluxes, although the spectrum below ~ 1 TeV be- 
comes softer with time. In Fig. [151 we present the evolution 
of the specific fiux at 5 GHz, and the bolometric fluxes be- 
tween 1-10 keV, 0.1 - 100 GeV, and > 100 GeV, up to 
10^ yr. In general, all the lightcurves show a smooth in- 
crease for the considered times. 

From the morphological point of view, the sources 
discussed here would appear extended at scales of ^ 
10' {d/ 3 kpc)~^ in radio, soft X-rays and high-energy 
gamma-rays. At hard X-rays and > 1 TeV, the emis- 
sion would be rather concentrated close to the jet re- 
verse shock and possibly the jet recollimation shocks, 
since the emitting electrons, with lifetimes ~ 1.3 x 
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Fig. 15. Evolution of the computed fluxes at 5 GHz (top, 
left), 1 - 10 keV (top, right), 0.1 - 100 GeV (bottom, left), 
and > 100 GeV (bottom right) up to 10^ yr, for cases 1 
(black/thick lines) and 2 (blue/thin lines). 

10^ {B/IQ-^G)-'^ (£;/10TeV)-^ yr, may not have time to 
fill the cocoon. 



5. Discussion 

We briefly discuss here, under the light of our results, the 
two HMMQs with very clear evidence of jet interaction with 
the ISM: SS 433 (case 1) and Cygnus X-1 (case 2). In the 
case of SS 433, the jet is strongly interacting with its su- 
pernova remnant (e.g. Zealey et al. 1980). This source is 
therefore a clear example of a jet/SNR interaction, but the 
luminosity of the jet, ~ 10^^ erg s~^, is much higher than 
the value adopted here. Nevertheless, there is evidence of 
non-thermal radio and X-ray emission in the SNR/jet in- 
teraction regions, which implies that particle acceleration is 
taking place there (e.g. Safi-Harb & Petre 1999), and possi- 
bly gamma-rays, although at present only upper-limits are 
available (see Bordas et al. 2010b and references therein). 
It is worth noting, that the interaction of an SS 433-like jet 
with its external medium has been detected even in other, 
nearby, galaxies (e.g. Pakull et al. 2010). 

Concerning Cygnus X-1, there is a thin emission region 
that seems to trace the bow shock formed by the jet in an 
inhomogeneous ISM (the counter-jet is not detected; see, 
e.g., Gallo et al. 2005). A SNR/jet interaction can indeed 
be discarded, because the SNR apparently was very weak, 
as reported by Mirabel & Rodrigues (2003). The struc- 
ture is elongated, and thus a shocked wind/jet origin seems 
unlikely from observations, although the wind energetics 
(~ 10^^ — 10^^ erg ~^) plus medium inhomogeneity might 
be a possibility, albeit remote. On the other hand, the de- 
projected distance of the bow shock to the system is about 
10 pc (Gallo et al. 2005), enough for the jet to have al- 
ready left the shocked wind/ISM region. As noted above, 
the non-detection of emission from the cocoon region may 
be related to jet disruption. In such a case, radiation may be 
still generated from particles accelerated via stochastic or 



shear acceleration or the jet recollimation shock(s), which 
is expected even if the jet disrupts at farther distances (see 
Bordas et al. 2010a for a similar situation in extragalac- 
tic sources). The bow shock in Cygnus X-1 does not show 
clear signatures of the system proper motion, which can be 
explained by the low peculiar velocity of the source with 
respect to its birth association, Cygnus 0B3 (Mirabel & 
Rodrigues 2003), plus a relatively young jet age (Gallo et 
al. 2005). 

Other high-mass systems may be good candidates to 
study the interaction of the jet with its surroundings. 
Cygnus X-3, for instance, has a powerful jet, and radia- 
tion has been detected in the surroundings of this source 
(Heindl et al. 2003; Sanchez-Sutil et al. 2008), although 
the non-thermal nature of this emission and its relation to 
Cygnus X-3 are still to be confirmed. Hints of extended X- 
ray emission have been also found in LS I +61 303 (Paredes 
et al. 2007; see however Rea et al. 2010), but the micro- 
quasar nature of this source is disputed (see, e.g., Bosch- 
Ramon & Khangulyan 2009 and references therein). 

The number of powerful HMMQs interacting with the 
environment may be of the order of 10 in our Galax}{3 (see, 
e.g., Perucho et al. 2010). The small number of sources, and 
their preferred location in the inner regions of the Galaxy 
(Bosch- Ramon et al. 2005), makes their detection difficult 
due to strong local/foreground absorption of radio and soft 
X-ray emission. Diffuse hard X-rays, clustered in the accel- 
eration regions (recollimation and reverse shocks, and pos- 
sibly the bow shock generated in the jet-ISM interaction), 
may be better probes to unveil young HMMQs. Very high- 
energy radiation could be detected from jet/medium in- 
teractions by forthcoming Cherenkov instrumentation (e.g. 
CTA), which in addition does not suffer of strong gamma- 
ray background contamination. It could be however that 
most HMMQs would be old and fast enough to have a well 
developed bow shock driven by the proper motion of the 
system. In such a case, the jet would be destroyed and 
mixed with material from the system motion-driven bow 
shock, and its radiation, if any, could be hard to detect. 

6. Summary 

In the present work, we explore different possibilities for the 
propagation of a HMMQ jet assuming different initial times 
for the jet activity to start with respect to the formation of 
the HMXB: shortly (40000 yr) after the SNR explosion of 
the progenitor of the compact object; at epochs (~ 10^ yr) 
at which only a spherical wind/ISM shock is present; and at 
the stages when the system peculiar motion bow-shapes the 
wind/ISM shock creating strong shocked wind convection 
backwards (> 10^ yr). For the two first cases, the jet needs 
about 10^ yr to reach the unperturbed ISM, a period much 
smaller than the potential duration of a HMMQ, which 
may be up to the lifetime of the remaining massive star, 
of ~ 10^ yr. For the third case, the jet gets disrupted and 
mixed with the shocked stellar wind, and advected away 
from the stellar bow shock, in just ^ 10^ yr. For the cases 
in which the jet reaches the normal ISM, its propagation 
can be described as in Bordas et al. (2009), although the im- 

^ We do not consider here the impact of microquasars as 
a class, including low-mass sources, on the ISM properties 
(see, e.g., Heinz et al. 2008), but focus on the specific case of 
jet /environment interaction in HMMQs. 
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pact of crossing the shocked wind/SNR eject a/ shocked ISM 
can lead to jet disruption, as in FRI jets (e.g. Perucho & 
Marti 2007). When the jet is embedded in the hot shocked 
wind/SNR ejecta, significant non-thermal emission can be 
generated in the cocoon region, mainly in radio and X-rays, 
but also gamma rays. Hard X-rays and very high-energy 
photons may be the best probes to unveil the presence of 
microquasar jets interacting with their environment. 
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